Localization of Spinwaves in the Quantum Hall Ferromagnet. 
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The quantum Hall ferromagnet at filling fraction v = 1 has some unusual properties due to the 
remarkable identity between the topological density of a spin distortion and the associated electrical 
charge density. We investigate the localization of spinwaves by coupling to a scalar disorder potential 
via their topological density. A low energy description of the system in terms of a non-linear sigma 
model of unitary supermatrices is derived. All states of this model are localized in two dimensions. 
A possible experimental signature of these effects in photo-luminescence is suggested. 



The ground-state of a two-dimensional electron 
gas(2DEG) at exact filling of a single Landau level (filling 
fraction v = 1) is strongly ferromagnetic. The properties 
of this quantum Hall ferromagnet(QHF) are profoundly 
affected by the topological nature of the quantized Hall 
plateau. In a quantum Hall state, there is a commensu- 
ration between the magnetic flux through the 2DEG and 
the electrical charge density. This leads to the identifica- 
tion of the topological charge density of a spin distortion 
with an associated electrical charge densitya (the Berry 
phase induced by a spin distortion may be reproduced by 
the Aharonov-Bohm phase induced by a fictitious mag- 
netic flux). As a consequence, the elementary excitations 
formed as the filling fraction is moved slightly away from 
v = 1 are electricallM^jind topologically charged objects 
known as SkyrmionsBH. In addition, spinwaves may cou- 
ple to a scalar disorder potential via their topological 
densityd. The geometrical structure of this disorder in- 
teraction is very similar to minimal coupling of quantum 
particles to a random flux and the spinwave system pro- 
vides a novel realization of this intently studied problcrnQ. 
In this paper, we investigate the localization of spinwaves 
in circumstances where the disorder potential is suffi- 
ciently weak that the ground-state remains ferromagnetic 
and coupling to spinwaves is the dominant effect. In the 
first section, we use diagrammatic techniques to demon- 
strate a decoupling of charge-density fluctuations from 
fluctuations in .exchange-energy. In the second, we use 
supersymmetryd to construct a low energy description of 
the system in terms of a non-linear sigma model of uni- 
tary super-matrices. All states of this model are localized 
in two dimensionsu. A possible experimental signature of 
these effects in photo-luminescence is suggested. 

Our starting point is a spinwave expansion about the 
continuum field theory of the QHF proposed by Sondhi et 
aoH. The effective action and electrical current density 
for small fluctuations, 1 = ia? 0)> about the ferromag- 
netic ground-state, n = (0,0, 1), are given by 



5* = 



d 2 xdt 



~f (~ ^ - p s V 2 - pgB) I ~ J (x)[/(x) 



J„ = i—e^ x djdxl. 
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tern. The first part of Eq. ([!]) is the usual Schrodinger ef- 
fective action for spinwaves in a continuum ferromagnet. 
We use the complex notation, I = l\ + Hi, I = li — Hi. 
p is the electron density (p = v/2nl B , where Ib is the 
magnetic length), p s is the spin stiffness and g is the 
Zeeman coupling, into which we have absorbed the elec- 
tron spin and the Bohr magneton for ease of notation. 
The second term Eq.(|l|) is due to the identity of charge 
and topological charge embodied in Eq.(||). It describes 
the interaction of spinwaves with a scalar disorder poten- 
tial, £7(x). The Coulomb self- interaction of the spinwave 
charge density has been neglected, since it is higher order 
in both derivatives and in the spinwave field. Through- 
out the bulk of the calculations presented in this paper, 
we will ignore the Zeeman term in Eq. (Jl|) . It is a simple 
matter to re-introduce it at the end of our calculations. 

The correlations in the disorder potential felt by the 
two-dimensional electron gas in GaAs heterostructures 
are conveniently modeled as followsiQ 



((<7 q C/ q ,)) = (2^(q + q')7- 



=-2|q|< 



(3) 



d is the width of the insulating spacer layer separating 
the electrons from the ionized donor impurities. 7' is a 
measure of the disorder strength and is related to the area 
density of donor impurities, rid; l' = (e v / nd/2e) 2 (^/87r) 2 . 
Our aim in this paper is to determine the universal fea- 
tures that arise due to the geometrical structure of the 
interaction. To simplify our explicit calculation, we as- 
sume a Gaussian, 5-function correlated distribution for 
the disorder potential: 



«0qt>q'» = (27r)V(q + q') 



(4) 



n = {h,h, yl — |1| 2 ) is an 0(3)-vector order parameter 
describing the local polarization of the quantum Hall sys- 



We return to the more realistic correlations of Eq.(j|) 

towards the end of the paper. 

Laid out in this way, the interaction of spinwaves with a 

disorder potential has many similarities to the problem of 
QYnon-interacting Schrodinger particles in a random fluxta. 

'To see this, one should integrate the interaction in Eq.(Q) 

by parts. The result looks like minimal coupling to a 
^ -vector potential, A; = — (eu/16Tr)cijdjU(x), aside from 

the absence of an |A| 2 term (which leads to broken gauge 

symmetry in the present case). 
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The self-energy, calculated in the self-consistent Born 
approximation, with correlations in the disorder poten- 
tial given by Eq.ffl), is 



I. DIFFUSONS AND COOPERONS 

The diffusive and weak localization effects in a disor- 
dered system are determined by the diffuson 



p+k, £i+(D 
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i j 



(5) 



In this diagram, full lines represent spinwave propoga- 
tors, cross-bars indicate spatial derivatives of these pro- 
pogators and the dotted lines represent disorder correla- 
tions. The thick line represents a full spinwave propoga- 



tor including the self-energy; G 



Gn 



E. The real 



part of the self-energy may be absorbed into a renormal- 
ization of the spinwave stiffness. The imaginary part is 
given by 

ImS fl (p, W )= 7 Im J ^0(p x q) 2 G*(q, W ). (6) 

To lowest order in the disorder strength, this has the 
solution 



ImS R (p,w) 
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The disorder averaged spinwave Green's function is then 



(G « (P ,„ )) = uw-f +4 & (f)iP 
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In writing down the hnal expression for the disorder av- 
eraged Green's function, we have removed a factor of 
(f + i'jpuj/Wp'^j and Taylor expanded the denominator 
in powers of jpto/Wp^.. Since this latter quantity is a 
small parameter in our perturbative expansion, this Tay- 
lor expansion is justified. 

The most important point to notice here is that the 
scattering time diverges as the frequency goes to zero. 
The spectral weight is concentrated in an energy range 
1/2t w of the bare pole at pu/2 — p s |q| 2 , such that 
1/2t w <C pui/2. This should be compared to the case 
of electrons scattering from a random scalar potentials, 
where the bare pole of the Green's function is near to 
the Fermi energy, Ep, and the scattering rate, 1/r, is 
constant. The validity of the perturbative expansion for 
electrons depends upon the smallness of the parameter 
1/(tEf). In the present case, the extra derivatives in 
the interaction vertex lead to the existence of a small 
expansion parameter without the existence of a Fermi 
surface. From a calculational point of view, pcu/2 plays 
the role of a frequency dependent Fermi surface and mo- 
mentum integrals may be carried out in the same way as 
one would in the case of electrons. 



D(p,q,k,fi,n + w) = 
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and Cooperon modes 



C(p,q,k,0,ft + w) 
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respectively. Notice that the interaction vertex is anti- 
symmetric under time reversal and, therefore, occurs 
with the opposite sign in the ladder summation for 
the Cooperon, compared with the Diffuson. The lad- 
der approximation to these diagrams, as indicated in 
Eqs.(|[fO|), may be summed to give two Dyson's equa- 
tions that may be represented by the following diagram: 




In fact, these equations are rather difficult to handle. The 
trouble is that the geometrical factors, involving cross- 
products of momenta, are not easy to factorize. One way 
to negotiate these difficulties is to write the interaction 
vertex as the sum of four vertices; 



i j 

-h+- 
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(12) 



The diffuson and Cooperon diagrams are decomposed 
similarly into the sum of diagrams; 
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D = l>,,r , + D xxxy + D xxyx + 13 other terms, 



(13) 



where the suffices label the derivatives on the external legs of the contributing diagrams. We wish to evaluate these 
diagrams in the limit of small momentum, k, and frequency, lo. Consider the momentum integral over some internal 
rung of a ladder diagram when k = 0; 



a b 



d 2 p 
(2irf 



PaPbPcPdG{p, 51 + U)G(P, O). 



(14) 



The angular integration gives zero if a =/= b = c = d (or a permutation of this). The remaining non-zero rungs are 



xx yy y x x y 

yy yy yx yx 



yy xx xy yx 

Only elements of D corresponding to these rungs are non-zero. In addition, rotational symmetry requires that 

Dxxxx — Dyyyyi D XX yy — ^yyxx ? 

D X yxy — Dy X y X: D X yyx — D yX xy- (1^0 

Using these symmetries in Eq.(0), with k = 0, the Dyson's equation decouples into two pairs of coupled equations, 



p,n+co 






(17) 



relating D yy yy to D xxyy and 



p, n+co 




x y x y 





(18) 
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relating D xyyx to D yxyx . The components, D xyxy , etc. 
of D, describe the propogation of charge-density fluc- 
tuations (the fields on the left-most legs occur in the 
combinations d\ldil and d<ild\l, the charge-density be- 
ing proportional to eijdildjl) and the components D xxxx , 
etc. describe the propogation of exchange-energy fluctu- 
ations (the fields on the external legs occur in the com- 
binations d\ld\l and 821821, the exchange-energy being 
proportional to dildil). The decoupling of the Dyson's 
equation, therefore, amounts to decoupling of charge and 
energy propogation. An alternative way of separating a 
Dyson's equation such as Eq.(pd|) is presented in Ref.[ |). 
Written in the above manner, the physical cignemodes of 
the scattering process are apparent. Using the ansatz 

D a bcd{p, q, k, ft, to) = (p + k) a p b (q + k) c q d D abcd {\i, w, ft) 

(19) 



Eq.(p7|) reduces to 

D xxyy — 7 + D xxyy (cos 2 9 sin 2 9) e + D yyyy (cos 4 9) g 
d 2 q 



X7 



(27T)' 



D 



yyyy 



Dxxyy (COS 



X7 



d 2 q 

(27T) 2 



|q| 4 G A (q,fi)G i? (q,^ + c) 



Dyyyy (COS 2 9 SUl I 



q\ 4 G A (q,fl)G H {q,fl + uj) (20) 



Evaluating the momentum integrals and performing the 
angular averages, one obtains 

DxXyy ( 1 - 1 ~ J = — + Dyyyy (\ + i2WT f2 ) 



Dyyyy I 1 i 



D r 



(1 + i2ur n ) . (21) 



To lowest order in to, these equations have the solution 



D xxyy (k = 0, to, ft) = D yyyy (k = 0, to, ft) = (1 2 , 



(22) 



D xxyy and D yyyy , therefore, describe massless diffusive 
modes. A similar calculation for D xyyx and Dy X y X pro- 
duces a result that is finite as to — > 0. These correspond 
to massive modes. Expanding the Dyson's equation to 
lowest order in k, one obtains the result 

D XX yy(\S., UJ, ft) = Dyyyy(k, bJ , ft) = D(\S., UJ , ft), (23) 

where 



D(k,oj,fl) = 2^{z^-+D \k\ 



(24) 



and Do = Psftrn is the classical diffusion coefficient for 
the spinwaves. In principle, one should worry that, when 



k ^ 0, the Dyson's equations no longer decouple into 
two pairs of equations.i-.This may lead to a shift in the 
diffusion constant, Z?oQi3. We find that, in this case, the 
result is unaffected to lowest order. Finally, summing all 
of the massless contributions to the diffuson, we find 

D{p, q, k, fl, ft + uj) = (p + k).p(q + k).q£>(k, w, ft), 

(25) 

which is the massless propogator for exchange energy 
fluctuations. The Cooperon diagrams may be evaluated 
similarly. They are finite in the limit k, uj — > and, 
therefore, massive. This is a consequence of the time re- 
versal anti-symmetry. There is no contribution to weak 
localization from these modes. 



II. SUPERSYMMETRIC SOLUTION. 

We now develop a low energy theory for the interaction 
of spinwaves with a weak disorder potential, using super- 
symmetric techniques. The main subtleties of the current 
problem are in the handling of the geometrical factors in 
the interaction. We will concentrate upon these features 
and refer the reader to the literatureo for details of the 
super-symmetry itself. As noted previously, this prob- 
lem is very similar to that of non-interacting particles in 
a random flux. An alternative derivation of the results 
^resented here may be made using the techniques of Ref . [ 

We wish to determine the disorder averages of dynam- 
ical quantities involving (G A G R ) and, therefore, intro- 
duce an eight component superfield, 

^(iWWxY), (26) 

in the usual wayi. The superscripts, A/R, label advanced 
and retarded sectors, the fields x are anti-commuting and 
the overbar indicates complex conjugation. The presence 
of the commuting and anti-commuting fields in ip removes 
the need to write the partition function explicitly in the 
denominator of correlation functions and allows the dis- 
order average to be performed immediately. The result- 
ing Lagrangian for the superfield, tp, is 



£ = 



-ir ( -p s V 2 - 



ip + 1 (tijdiipdjipy 



d 2 r, 



(27) 



where e = (fl+uj/2)l + (ui/2 + i5)A and A is the diagonal 
supermatrix with elements ±1 in the advanced/retarded 
sectors. As is usual in the derivation of the supersym- 
metric sigma model, we assume that to <c ft. Disorder 
averaging may only induce cross correlations between the 
Green's functions if the difference in frequencies is less 
than the scattering rate puj/2 < 1/2tq,1/2tq +u1 <C pft/2, 
therefore this approximation is justified. Next, we decou- 
ple the quartic interaction, introduced by the disorder av- 
erage, with a supermatrix field, Qij, where i,j £ {x,y} 
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and each element of the 2x2 matrix Qij is an 8 x 8 
supermatrix; 



cxp (-A nt [ip, Qij]) DQi 



S = 



P 




P ( P^ 



d 2 rStr [Q 

2 ,, 



dp f d q 




The supertrace, Str, is as defined in Ref. [ ^| and summa- 
tion over repeated spatial indices is implied. Integrating 
out the superfield, ip, we obtain the following free energy 
functional for Qif. 



2p s J J (2tt) 2 7 (2tt) 2 
x Str [e ki pj (p + q)kG(p)5Q i:j (q) 

x e n i(p + q) m p n G(p + q)8Qi m (-q)] . (34) 

p- integral is only finite for the products G A G R , 
which are generated by the components of Q off-diagonal 
in the advanced and retarded sectors. Henceforth, we use 
Q to indicate supermatrices with only these off-diagonal 
components non-zero. Carrying out these integrals and 
substituting Q s , Q c \ and Q C 2 from Eq. (|33"|) , we find the 
following effective action for SQ S : 



F[Qi 



Str 



(lnG- 1 ) + -^— Q 
2 64tq,o s 



d 2 , 



(29) 



where G(r,r',Q) is the supermatrix Green's function of 
the field ip and satisfies the equation 



S = 



where 



1 



d 2 q 



32p s J (2tt 



-ipw/2 + D \q\ 2 ) 

xStr[5Q s (q)6Q s (-q)}, (35) 



pCl' 
2p~s 

xG{t,t',Q) = itf(r-r'). 



£fci [Qijdjdk + d k Qijdj]J 
(30) 

The saddle point equation for is 

f cPp 



Do = PsTntt. 



(36) 
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;£ktPkPjG{v : Q l:j ). (31) 



This is precisely the self-consistency equation that was 
solved previously to find ImS. The solution is 



eijVAV, 



(32) 



The corresponding actions for SQ C \ and SQ C 2 contain 
massive propogators. It is important that one should not 
simply ignore these massive modes. They may lead to a 
renorrualization of the diffusion constant for the massless 
modesBu. Here, we find that this is not the case. The 
cross-terms between 5Q S and SQ c i/SQ C 2 are proportional 
to |q| 2 . Integrating out the massive modes induces a |q| 4 - 
term in the SQ S propogator, which we neglect at small 
momentum. Keeping only fluctuations of Q over the sad- 
dle point manifold, Eq. ([35|) determines the coefficients of 
the terms in the sigma model; 



where V is an arbitrary, unitary supermatrix such that F\Q 1 = - Str / Dq\WQ 



VV = 1. The diagonal terms of Qij are zero at the sad 
die point and, therefore, correspond to massive modes. 
The off-diagonal components, however, sit in a Mexi- 
can hat potential and have massive longitudinal fluctu- 
ations and massless transverse fluctuations. At the sad- 
dle point, Qij ~ ekidkijOjip, and, therefore, the diago- 
nal elements of Qij describe charge density fluctuations, 
Qu ~ ^ijdiipdjip and the off-diagonal elements describe 
exchange energy fluctuations, tijQij ~ diipdiip. This is 
the same decomposition into massive and massless modes 
as was found diagrammatically in the preceding section. 
Let us define 

Qij = fijQs + diag(Q c i, <2 c2 ) 



128 Ps 



2 -2zf^ 



V 2 



AQ S 



d 2 r 



(37) 



Do is the classical diffusion constant given by Eq.(p6[). 
It may be written in the form Do — prnVQ/2, where 
vl = (dE/dp) 2 | 

E=n/2 ~ 2Ps^/p- This energy func- 
tional describes the diffusive propogation of exchange en- 
ergy deasity fluctuations. All states are localized in this 
modelM, with a localization length 



£ = vnm exp 



n 2 D 2 



64p 2 



rr 3/2 exp [ir 



(38) 



Qs 



i 



i Qi 



Qci = Qu 

Qc2 = Q22 



(33) 



At the saddle point Q s = VAV and Q c = 0. Expanding 
the free energy functional, Eq.(p9|), to quadratic order in 
fluctuations of Qij about the saddle point, we obtain 



which is divergent in the £1, q = limit. This divergence 
is a natural consequence of the existence of an SU(2) 
global symmetry (broken to U(l) with the inclusion of 
the Zeeman coupling). Inclusion of the Zeeman coupling 
modifies the parameters of the sigma model such that f2 
is replaced by Q—gB. There are no spinwave states below 
the Zeeman gap and the localization length diverges as 
the frequency approaches this, q = 0, limit. Realistic 
correlations in the disorder potential, given by Eq. ([?]), 
lead to a further modification of the parameters in the 
sigma model. The scattering time is given bycl 
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which again diverges as £1 — > such that 1/tq <C pO/2. 
Subsequent calculations follow through as before, with 
some additional complications. Finally, one obtains the 
same supersymmetric sigma model, Eq.(p7[), with 

D' = pSoCl = (40) 
7 

The diffusion constant is now independent of f2 and the 
localization length no longer has an exponential depen- 
dence upon frequency, but a power law dependence; 

'1/2 



(41) 



The numerical prefactor in this expression has an ex- 
ponential dependence upon the disorder strength. Esti- 
mates assuming a donor density, n c i, of the same order as 
the electron density, p, gives prefactors upwards of 10 3 . 
Experimental determination of the diffusion coefficient 
and localization length may be possible using space/time- 
resolved photo-luminescence. In the QHF, magneto- 
excitons and spinwaves are identical! If, in addition, 
the sub-band wavefunction of valence holes is centered 
in the same position as the wavefunction of electrons in 
the two-dimensional electron gas, then escitons have the 
same Hamiltonian as magneto-excitonstj. Localization 
and diffusion of excitons in other systems have been mea- 
sured using photo-emission spectroscopy^. The same 
techniques may be applicable here. 

Above a threshold disorder strength, the groundstate 
ceases to be ferromagnetic and contains a-xandom distri- 
bution of Skyrmions and anti-Skyrmionsalij. In certain 
circumstances, the main effect of this topologically non- 
trivial groundstate upon the spinwave action iSrto intro- 
duce an effective magnetic field and mass terir£3. Both 
the magnetic flux and the mass are proportional to the 
topological density of the background. This is yet an- 
other variant of the random flux problem. In this case, 
the effective random flux may be very large and localiza- 
tion may have further experimental consequences!!! 

In summary, we have considered the effect of weak dis- 
order upon the spinwave dynamics in the QHF. We find 
a spinwave scattering time that diverges as frequency 
approaches zero, thus ensuring the validity of perturba- 
tive and supersymmetric techniques despite the lack of 
a Fermi-surface. Charge and exchange-energy fluctua- 
tions arc found to decouple, the former being massive 
and the latter diffusive. The low energy dynamics of 
this system are described by a non-linear sigma model of 
unitary supermatrices. All states in this model are local- 
ized. A possible signature of these effects may be found 
in space/time-resolved photoemission spectroscopy. 

Thanks to Nigel Cooper, Michael Fogler, Duncan Hal- 
dane, Matt Hastings, Shivaji Sondhi and Alexei Tsvelik 
for comments and suggestions. 



t Present Address: T.C.M., Cavendish Laboratory, Mading- 
ley Rd., Cambridge, CB3 OHE, U.K. 

1 S. L. Sondhi, A. Karlhede, S. A. Kivelson and E. H. Rezayi, 
Phys. Rev. B47, 16419 (1993). 

2 S. E. Barrett, G. Dabbagh, L. N. Pfeiffer, K. W. West and 
R. Tycko, Phys. Rev. Lett. 74, 5112 (1995); R. Tycko, S. 
E. Barrett, G. Dabbagh, L. N. Pfeiffer and K. W. West, 
Science 268, 1460 (1995); E. H. Aifer, B. B. Goldberg, D. 
A. Broido, Phys. Rev. Lett. 76, 680 (1996); A. Schmeller, 
J. P. Eisenstein, L. N. Pfeiffer and K. W. West, Phys. Rev. 
Lett. 75, 4290 (1995). 

3 A. G. Green, Phys. Rev. B 57, R9373 (1998). 

4 A. G. Aronov, A. D. Mirlin and P. Woffle, Phys. Rev. B49 
16 609 (1994). 

K. Efetov, Super symmetry in disorder and chaos, C.U.P. 
(1997). 

6 K. B. Efetov, A. I. Larkin, and D. E. Khmelnitskii, Zh. 
Eksp. Teor. Fiz. 79, 1120 (1980) [Sov. Phys. JETP 52, 568 
(1980)]; L. Schofer and F. J. Wenger, Z. Phys. B 38, 113 
(1980); K. B. Efetov and A. I. Larkin, Zh. Eksp. Teor. Fiz. 
85, 764 (1983) [Sov. Phys. JETP 58, 444 (1983)]. 

7 A. L. Efros, Solid State Com 70, 253 (1989). 

8 P. Woffle and R. N. Bhatt, Phys. Rev. B 30, R3452 (1984). 

9 C. Kallin and B. I. Halperin, Phys. Rev. B 30, 5655 (1984). 

10 N. R. Cooper and D. B. Chklovshii, Phys. Rev. B 55, 2436 
(1997) 

11 H. F. Hess, E. Betzig, T. D. Harris, L. N. Pfeiffer and K. 
W. West, Science 264, 1740 (1994); Y. Takahashi, S. S. 
Kano, K. Muraki, S. Fuktsu, Y. Shiraki and R. Ito, Appl. 
Phys. Lett. 64, 1845 (1994); N.-H. Ge, C. M. Wong, R. L. 
Lingle Jr., J. D. McNeill, K. J. Gaffney and C. B. Harris, 
Science 279, 202 (1998). 

12 A. M. Polyakov, Gauge Fields and strings (Harwood Aca- 
demic, New York 1987). 

13 A. J. Nederveen and Yu. V. Nazarov, |cond-mat/9808339 

14 In preparation. 



G 



